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Low abundances of heavy elements in the solar outer layers: 
comparisons of solar models with helioseismic inversions 
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ABSTRACT 

Context. Recent solar photospheric abundance analyses have led to a significant reduction of the metal abundances compared to the 
previous determinations. {Z/X)s is 30% smaller, down to 0.0165. The solar models computed with standard opacities and diffusion 
processes using these new abundances give poor agreement with helioseismic inversions for the sound-speed profile, the surface 
helium abundance, and the convective zone depth. 

Aims. We attempt to obtain a good agreement between helioseismic inversions and solar models that present the "old mixture" in 
the interior and new chemical composition in the convective zone. To reach this result, we assume an undermetallic accretion at the 
beginning of the main sequence. 

Methods. We compute solar models with the Toulouse-Geneva Evolution Code, in which we simulate an undermetallic accretion in 
the early stages of the main sequence, in order to obtain the new mixture in the outer convective zone. We compare the sound-speed 
profile, the convective zone depth, and the surface helium abundance with those deduced from helioseismology. 
Results. The model with accretion but without any mixing process inside is in better agreement with helioseismology than the solar 
model with the new abundances throughout. There is, however, a spike under the convective zone which reaches 3.4%. Furthermore, 
the convective zone depth and the surface helium abundance are too low. Introducing overshooting below the convective zone allows 
us to recover the good convective zone radius and the addition of rotation-induced mixing and tachocline allows us to reconcile 
the surface helium abundance. But in any case the agreement of the sound-speed profile with helioseismic inference is worse than 
obtained with the old abundances. 
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c3 1. Introduction 



The standard solar models that include the metal abundances 
as given by Grevesse & Noels (1993) (hereafter GN93) and 
. the current physics updates (new opacity tables and equations 
' of state frorn OPA L, and element diffusion and settling, e.g., 
' iRichard et al.ll2004l) are in very good agreement with helioseis- 
mology. The square of the sound speed derived from this stan- 
dard solar model agree s with the values obtained from heli oseis- 
mology within 0.5% dChristensen-Dalsgaard et aT] 1996h. The 
recent solar photospheric abundances derived by 'Asplu nd et alJ 
(l2005h (hereafter Asp05) revise the CNO and Ne abundances 
downwards by more than 30% compared to the earlier determi- 
nations. The solar photospheric Z/X is reduced to 0.0165 + 10% 
and Z to 0.0122, whereas the "old" mixture gave a higher Z/X of 
0.0245. New solar standard models including the Asp05 mixture 
lead to a large disagreement with the helioseismically inferred 
sound-speed profile, convective zone depth, and surface helium 
abundance Ys . 

Several research groups have explored different ways to try 
restoring agreement between solar models and helioseismology. 
[Bahcall et al. ( 2005a) showed that an increase of 21% of the 
OR4L opacities could reconcile the new solar abundances with 
the helioseismic measurements of the depth of the convective 



zone. The uncertainties in the ORA^L96 opacity tables however 
are estimated to be no more than 4%. The uncertainties on the 
element diffusion velocities below the convective zone are larger 
and can reach 35% for Fe and O, the two main c ontributors to 
the opacity in this region dMontalban et al.l [20061) . An increase 
of the diffusion coefficient by a factor of 1.5 to 2 combined 
with a small increase of the opacity at the bottom of the convec- 
tive zone provides better agreement with helioseismolog y for the 
depth of the convective zon e and the profile tMontalban et al.l 
12004 iBasu&Antial 12004 . but in this case Ys has to be de- 
creased below the value de termined by the inversion of solar 
oscillations. The attempt bv iGuzik et alj (|2005.) to increase the 
diffusion coefficients for helium and the other elements sepa- 
rately to fit the Asp05 mixture in the outer layers and recover 
the GN93 mixture at t he bottom of the c o nvect ive zone leads 
to the same problem. Turck-Chieze et all (|2004|) showed that 
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the ne w oxygen measurement of the Sun from lAsplund et aT] 
(I2005h was in accordance with the observations of extragalac- 
tic HII regions, Magellanic clouds, and other clusters. They up- 
dated their predicted neutrino fluxes using the new composition 
and they obta i ned a goo d agreement with the detected values. 
lAntia & Basul (l2005h and lBahcall et al.l (l2005bh computed mod- 
els with an arbitrarily increased Ne/O ratio and showed that they 
could be in agreement with helioseismology. At the same time, 
Drake & Testa (20(j5h assessed a n higher Ne /O ratio for the Sun 
than the one adopted in|A splund et alj ( 2005h and consisten t with 
the results of lAntia & Basu (.20051) and bahcall et alj ( l2005bi) . by 
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observing several cool nearby stars in X-ray spectra. However, 
new measurements by analyzing the solar active region spectra 
(Schmelz et al. 2005) and the extreme ultravio let emission lines 
of NelV-VI and OIII-V ions in the solar corona (lYoung|2005l) . are 
in agreement with the Ne/O ratio value of the new abundances. 

Here we compute solar models with the assumption that un- 
dermetallic accretion occurred in the very early stages of the so- 
lar life on th e main sequence. A similar idea was also proposed in 
iGuzik et al.l ( I2005D without computing any model. We suppose 
that, during the process of planetary formation, the young Sun 
was surrounded by metal-poor matter while metals were concen- 
trated in planetesimals, and that an important part of this metal- 
poor matter fell onto it. We do not justify this idea here by any 
computations of accretion processes: our aim was only to com- 
pute "exotic" solar models that could be able to reconcile the 
new solar abundances with helioseismology. We simulated this 
undermetallic accretion to obtain the Asp05 mixture in the outer 
convective zone with an increased abundance inside. In this way, 
we can obtain a metallic difference between the outer layer and 
the inside larger than possible from diffusion processes only. In 
these models, we attempted to recover the sound-speed profile, 
the depth of the convective zone, and the surface helium abun- 
dance as given by helioseismology. 

In Sect. 2, we discuss the new Asp05 abundances and their 
influence on the sound-speed profile. In Sect. 3, the models with 
accretion are presented. The effects of overshooting below the 
convective zone and mixing induced by the tachocline and by 
rotation are respectively discussed in Sects. 4 and 5. The conclu- 
sions are given in Sect. 6. 
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Fig. 1. Comparison between the profile of sound velocity of the 
different models (SI : solid line, S2 : dotted line, S3 : dashed 
line, S4 : dot-dashed line, S5 : three dot-dashed line, S6 : long 
dashed line) and the one deduced from the helioseismology 



2. Solar models with the new Asp05 mixture 

The models are computed with the Toulouse - Genev a 
Stellar Evolution Code (TGEC) (iRichard et alj Il996h . 
We use the equa tion of state table s of OPAL2001 (up- 
dated versio n of Rogers et al.l 19961) . the opacities tables 
of OPAL dRogers& IglesiasI Il995h completed with the 
[Alexander & Fergusson' (1994') low temperature opaciti es, and 
the N ACRE compilation of nuclear reaction rates ( Angulo et alj 
11999 ) with the Bahcall screening routine. The convection treat- 
ment is based on the mixing-length theory dBohm-Vi tense* 1 958l) 
and t he diffusion coefficients are computed as in iPaq uette et al. 
([19861) . 



Model SI is a standard model, computed with the "old" 
GN93 chemical composition. No mixing other than convec- 
tion is introduced. The calibration procedure includes two free 
parameters: the initial mass fraction of helium Yq and the 
mixing-length parameter a. These parameters are adjusted to 
obtain the solar luminosity and r adius at the so l ar age : Lq- 
3.8515 + 0.0055 x 10^^ erg.s"' (lGuenthereta lJ ll992h R^^ 
6.9575 ± 0.0024 x 10^° cm (mean obse rved value : lAllenlll976t 
iBrown & Christensen-Dalsgaardl[T998 ^. and = 4.6 Gyrs. The 
characteristics of model SI are presented in Table[T] 

This model is in good agreement with helioseismology. The 
comparison of the sound velocity computed in this model with 
the one deduced from helioseismology by iBasu et al. (1997) 
is presented as the solid line in Fig. [1] The relative difference 
between model S 1 and the seismic inversions is always smaller 
than 0.5%. The depth of the convective zone is also consistent 
with helioseismology: we find rgz/Re - 0.712 compared to 
r,-JRa = 0.713 ± 0.001 (Bas u & Anti"alll997h . The helium 
abundance is 0.240, slightly smaller than the value obtained 
from seismic inversions, which is typically 0.249 + 0.002 



(lBasu& Antialll997 ^. This is due to the fact that model SI does 
not include any mixing process which could slightly slow down 
hehum settling. 

Model S2 has been computed using the Asp05 mixture 
throughout. Its characteristics are presented in Table [1] and the 
sound velocity in Fig.[T] The discrepancy between the model and 
the seismically derived sound velocity is much larger than for 
model SI: it goes up to 3%. The initial Z/X is lowered by 33%, 
compared to model SI: the new initial value is 0.0183 instead 
of 0.0272. More precisely, the CNO and Ne mass fractions have 
decreased by 36%. The quantitative analysis of iTurcotte et al.l 
(Il998l) for the GN93 mixture shows that the main contributors 
to the radiative opacity below the solar convective zone are O, 
Fe, H, Ne, and, with smaller contributions, N and C. Hence, the 
reduction of the metal mass fraction and the decrease of O and 
Ne abundances induce important changes in the opacity related 
to the Asp05 mixture: the opacity at the bottom of the convec- 
tive zone decreases by 25% compared to the GN93 case. As the 
radiative temperature gradient V^rf is proportional to the opac- 
ity K, the new convective zone is shallower (JczIRq - 0.729). 
Furthermore, the helium abundance in the convective zone is 
lower than the helioseismic value (I5 = 0.223). 



3. A new undermetallic model 

The problem of new solar models seems very difficult to solve, 
as the new determinations of solar abundances, using new tech- 
niques of 3D-modelization, are incompatible with the seismic 
studies of the Sun. On the other hand, the old GN93 solar mix- 
ture gave a very good consistency with the helioseismic analysis, 
within 0.5%. We decided here to compute new models using the 
GN93 mixture in the solar interior, to try to remain consistent 
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Calibration 
parameters 


External 
parameters 


Convective 
zone 




a 




L(10" erg.s"') 


R (10"'cm) 




Z/X 




Model SI 
(GN93) 


1.8104 


0.2716 


3.8517 


6.9575 


0.240 


0.0244 


0.713 


Model S2 
(Asp05) 


1.6556 


0.2562 


3.8511 


6.9576 


0.223 


0.0164 


0.730 


Model S3 
(GN93 
+ accr) 


1.6338 


0.2720 


3.8514 


6.9586 


0.240 


0.0165 


0.732 


Model S4 

(GN93 
+accr +ov) 


1.6344 


0.2721 


3.8519 


6.9552 


0.243 


0.0166 


0.712 


Model S5 
(GN93 
+accr +mix) 


1.5101 


0.2691 


3.8514 


6.9592 


0.249 


0.0167 


0.751 


Model S6 
(GN93 
+accr +mix +ov) 


1.5495 


0.2705 


3.8515 


6.9593 


0.249 


0.0165 


0.712 



with helioseismology, but with the Asp05 chemical composition 
inside the convective zone. This situation cannot be justified in 
terms of element diffusion only, as the precisio n on the diff'usion 
veloc i ties do es not allow such a large settling dMontalban et al.l 
12004 12006'). For this reason, we assumed accretion of metal- 
poor gas at the beginning of the main sequence. 

We simulated this accretion in a simplified way, by instanta- 
neously decreasing the metal abundances in the convective zone 
and increasing hydrogen and helium accordingly, in one of our 
first main sequence models (aged 74 Myrs). If tz is the factor 
(smaller than one) by which the abundances of metals are mul- 
tiplied, the hydrogen and helium abundances are multiplied by 
TxY, SO that : 



TXY - 



l- (l-X-Y)*Tz 
X+Y 



(1) 



In this simplified way, all metals are decreased with the same 
factor, so that the mixtu re remains unchanged , while the overall 
ratio (Z/X)s reaches the lAsplund et alj (l2005h value. 

Model S3, whose characteristics are also presented in Table 
[1] has been computed with tz = 0.5. The comparison of the 
sound velocities with the seismic inversions is presented in Fig. 
[U The discrepancy is smaller than for model S2, but a spike 
remains below the convective zone, which reaches 3.4%. This 
spike can be due to the fact that the convective zone is too shal- 
low (see Table [1) and also to the metal gradient induced by 
accretion and diff'usion below the convective zone (see Fig. |2]i. 
Meanwhile, the helium abundance in the convective zone is still 
too low. Two kinds of improvements are possible at this stage: 
introducing an overshooting zone below the convective zone so 
that the combined mixing reaches the depth deduced by helio- 
seismology (rcz/Re - 0.713), and adding a rotation-induced 
mixing below the convective zone to smooth down the metal gra- 
dient. In the following, we present three new models, one with 
overshooting (S4), one with rotation-induced mixing (S5), and 
one with both processes (S6). 

4. Convective overshooting 

The introduction of overshooting below the convective zone al- 
lows us to simulate an overall mixing zone down to the depth 
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Fig. 2. Metal profile below the convective zone of models S3 
(dashed line), S4 (dot-dashed line), S5 (three dot-dashed line), 
and S6 (long dashed line) 



deduced from helioseismology. This result is obtained with an 
overshooting parameter of 0.25 Hp, where Hp is the pressure 
scale height. Due to overshooting, a larger value of the accre- 
tion parameter is needed to obtain the Asp05 abundances in the 
convective zone at the solar age: tz = 0.58. The characteristics 
of the resulting model S4 are presented in Table [T] We can see 
that the bottom of the convective-plus-overshooting zone is the 
same as in model SI, namely r/R^ = 0.712. The comparison 
of the sound velocities with helioseismology is presented in Fig. 
[T]and the metal gradient in Fig.|2] The amplitude of the spike 
under the convective zone is smaller than without overshooting. 
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but it still reaches 2.1 
too low. 



The surface helium abundance is still too much compared to the observed value. 



5. Rotation-induced mixing 

Mixing below the convective zone was not introduced in the pre- 
vious models. As a consequence, lithium was not depleted as 
observed in the Sun. Solar spectroscopic observations compared 
to abundance determinations in meteorites show that the solar 
lithium has been depleted by a factor of 160 compared to its 
initial abundance an d that berylliuni has not been destroyed by 
more than 0.09 dex (lAsplundll2004l) : these features have to be 
taken into account in consistent solar models. Rotation-induced 
mixing can deplete lithium as observed and it also smooths down 
the metal gradient induced by accretion. We computed model 
S5 (see Table [1]) without overshooting, but with the rotation- 
induced mixing calibrated to have the suitable lithium destruc- 
tion at the solar age in qu estion. We used the prescription by 
iTheado & Vauclairl (l2003h . with a diffusion coefficient D,urb, 
such that: 



dc 1 d 



dt 
where 



dc 



iitrb 



dr 



c,.+ 



30C,, 



dr 



r\Ur\ = auirbAUri 



(2) 



(3) 



Here \Ur\ represents the vertical component of the meridional 
circulation velocity, and C/, and Cy are characteristic parameters 
for the horizontal and vertical diffusion coefficients: 



Dh = ChAUA 



(4) 
(5) 



Ch and Oturb are the free parameters used for the calibration that 
determines the efficiency of th e turbulent motions. We use the 
values given bv lTheadol d2002l) : C/, = 50000 and an,rb = 6. 

We also include in these computations a tachocline cali- 
brated by two coefficients: the diffusion coefficient at the base 
of the convective zone Dbcz and the thickness of the tachocline 
A. The diffusion coefficient inside the tachocline is given by: 



D,acho = £>teexp ln2 



(6) 



where is the radius at the base of the convective zone. 
These parameters are adjusted to reproduce the observed solar 
lithium depletion at the solar age: Dbcz - 8-0 x 10^ cm^.s"' and 
A - 0.048 X 10'° cm. As mixing reduces element diffusion, a 
stronger accretion is needed. In model S5, the accretion factor 
was taken as : tz = 0.37 and it was introduced in seven time 
steps. The corresponding metal profile in the present Sun is 
shown in Fig. |2l The comparisons of the sound velocity in 
model S5, compared to the seismic value, is presented in Fig.[T] 

The agreement is worse than for previous models. The spike 
under the convective zone reaches 5.7% and the difference 
with the seismic value remains large far below the convective 
zone, down to rjRc, - 0.3. This is due to the spread of the 
metal gradient in the interior and the consecutive changes 
in the opacity. Here the convective zone is still too shallow 
{vczIRq - 0.751), but the helium abundance is correct in the 
convective zone (Ys - 0.249) as well as the lithium value 
(Li/Lio - 1/167). Beryllium is depleted by 0.4 dex, which is 



Model S6 has been computed with both overshooting and 
rotation-induced mixing. The accretion and mixing parameters 
were modified compared to previous models : tz - 0.48 intro- 
duced in three time steps, Ci, — 65000, a„irb — 1, Dbcz ^2.5x10^ 
cm^.s ', and A - 0.048x10'" cm. The sound-speed profile com- 
pared with helioseismology is presented in Fig.[T]and the metal 
gradient in Fig.|2] The metal gradient is steeper in model S6 than 
in model S5 because, thanks to the overshooting, the rotation- 
induced mixing which reproduces the lithium underabundance 
has to be less efficient. Here beryllium is only depleted by 0.12 
dex. The overshooting allows us to reach the right depth of the 
convective zone, and the mixing gives the right surface helium 
abundance, but the sound-speed profile is still in large disagree- 
ment with helioseismology. 

6. Conclusion 

The incompatibility between the new solar photospheric abun- 
dances of As plun d et al. (2005) and helioseismology have al- 
ready been studied by many authors (B ahcall et al. 2005^ 
Montalban et"!!] |2004 120061 iBasu & Antia .2004; iGuzik etal] 
20051) . No satisfying solution has been found. The best agree- 
ment is obtained by increasing the opacity below the con- 
vective zone and the diffusion rates with suitable factors 
(IMontalban et al.l2006l) . But the required changes are larger than 
the accepted uncertainties, and the models still have too low sur- 
face helium abundances. In the present paper, we have computed 
solar models with the GN93 mixture in the solar interior and the 
Asp05 abundances in the external convective zone. These mod- 
els could be justified by accretion of metal-poor gas in the early 
stage of the main-sequence. 

Models in which such an accretion is introduced present a 
better agreement with helioseismology than models using the 
Asp05 mixture throughout, but there is still a large discrepancy 
below the convective zone in the sound velocity, and the depth of 
the convective zone as well as the surface helium abundance are 
incorrect. The introduction of overshooting below the convective 
zone allows us to reach a mixing zone as deep as that deduced 
from helioseismology. Adding rotation-induced mixing allows 
us to obtain the right surface helium abundance, and also the 
right lithium depletion. While the introduction of overshooting 
does not change the sound-speed profile in a significant way, the 
smoothing and spreading of the metal gradient due to rotation- 
induced mixing increases the difference with the seismic sound 
velocity profile in the interior We tried our best to reconcile 
the new solar abundances with helioseismology, but the results 
are not encouraging: it does not seem possible to obtain a solar 
model in good agreement with the seismic inversions, contrary 
to models computed with the old abundances. 
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